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Abstract. The millisecond pulsars discovered in globular 
clusters arc diagnostic tools for studying the dynamics of 
the clusters, the intracluster ionized gas and the evolution 
of the binaries embedded in the cluster. In the first two 
years of a search at 1.4 GHz in progress at the Parkes radio 
telescope 12 new millisecond pulsars have been discovered 
in 6 globular clusters in which no pulsars were previously 
known. The scientific fall-back of these discoveries is here 
discussed for the case of two globular clusters: NGC 6397 
and NGC 6752. 



1. Introduction 

Millisecond pulsars (MSPs) are old neutron stars recycled 
through transfer of matter and angular momentum from a 
mass losing companion in a binary system (e.g. Alpar et al. 
1982; Smarr & Blandford 1976; Bhattacharya & van den 
Heuvel 1991; Kulkarni & Anderson 1996). They are mas- 
sive and almost pointlike objects, appearing very good test 
masses for probing gravitational effects. Most of them are 
also extremely stable clocks, allowing very accurate mea- 
surements of their rotational parameters and position in 
the sky. Finally they are the endpoints of a complex stellar 
evolution in binary systems, involving not yet completely 
understood mechanisms like common envelope, accretion 
disk instabilities, magnetic braking, irradiation, evapora- 
tion. 

Besides evolution of a primordial binary, in a Globular 
Cluster (GC) another formation channel for the MSPs is 
available: in fact exchange interactions in the ultra-dense 
stellar environment of the cluster core can sustain the for- 
mation of various kinds of binaries suitable for recycling 
neutron stars (Davies & Hansen 1998; Rasio, Pfahl & Rap- 
paport 2001). As a consequence, despite the large differ- 
ence in total mass between the disk of the Galaxy and the 
Globular Cluster system, about 50% of the entire MSP 
population has been found in the latter. 

In particular millisecond pulsars found in globular clus- 
ter have proven to be valuable tools for studying: 

— the potential well of the globular cluster (e.g. Phinncy 
1992; Camilo et al. 2000); 



— the dynamical interactions in the globular cluster core 
(e.g. Phinncy & Sigurdsson 1991); 

— the gas content in the globular cluster (e.g. Freire et 
al. 2001) 

— the neutron star retention in the globular cluster (e.g. 
Rappaport et al. 2001); 

— the binary evolution in a crowded environment (e.g. 
Rasio et al. 2000). 

Unfortunately, the MSPs in GC are sources difficult 
to be detected. The main reason for their elusiveness is 
that GC-MSPs are often distant pulsars in close binary 
systems. Their large distances (i) make their fluxes typi- 
cally very small and (ii) strongly distort their signals due 
to the dispersive effects of propagation through the in- 
terstellar medium. Their inclusion in tight binaries (Hi) 
causes Doppler-shift changes of the apparent spin period 
and sometimes (iv) makes the radio signal periodically ob- 
scured by eclipses. 

As a consequence, except for the relevant case of 47 Tu- 
canae (Camilo et al. 2000), the discovery rate of these ob- 
jects strongly declined in the second half of the 1990s: in 
the 7 years from 1987 (when the first GC-MSP, B1821-24, 
was discovered in M28 at Jodrell Bank: Lyne et al. 1987) 
to 1994 (Biggs et al. 1994) 32 GC-MSPs entered the cata- 
log, whereas no new source was published in the following 

5 years. 

Since two years a new deep search has been running 
at the Parkes radio telescope; it involves people from the 
Cagliari Astronomical Observatory (N. DAmico), from 
the Australia Telescope National Facility (R. N. Manch- 
ester and J. Sarkissian), from the Jodrell Bank Observa- 
tory (A. G. Lyne), from the Columbia Astrophysics Lab- 
oratory (F. Camilo), and the author of this review. As a 
result, to date a dozen new pulsars have been discovered in 

6 globular clusters, none of which had previously known 
pulsars associated. More recently, 3 further millisecond 
pulsars have been detected in one of these clusters (namely 
NGC 6266) by Jacoby et al. (2002), exploiting the sensi- 
tivity of the new Green Bank Radiotelescope. So the total 
number of known millisecond pulsars in globular cluster 
is now 57. 
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Fig. 1. Observed signal intensity at 1.4 GHz as a function 
of orbital phase and pulsar phase for four long observa- 
tions away from the nominal eclipse region. The data are 
processed in contiguous integrations of 120 s duration, (a) 
~ 10 hr observation starting on 2000 December 3 at 22:23 
UT; (b) ^ 4 hr observation starting on 2001 March 1 at 
22:32 UT; (c) ~ 7 hr observation starting on 2001 March 
5 at 18:31 UT; (d) ~ 11 hr observation starting on 2001 
March 12 at 15:02 UT (from D'Amico et al. 2001c). 



In the following I will focus on the globular clusters 
NGC 6397 and NGC 6752, summarizing the main scien- 
tific results obtained so far. sometimes combining radio 
data with optical and X-ray observations. What I will re- 
view here have been originally reported in seven differ- 
ent papers: D'Amico et al. (2001b) and Possenti et al. 
(2001) described the new hardware and software systems 
adopted in the Parkes search for GC-MSPs, the current 
status of the experiment and its perspectives. D'Amico et 
al. (2001a) announced the discovery of the first 4 millisec- 
ond pulsars detected during the Parkes search for GC- 
MSPs. D'Amico et al. (2001c) and Ferraro et al. (2001) 
presented the timing solution and the optical observations 
for the binary millisecond pulsar in NGC 6397. D'Amico 
et al. (2002) and Colpi, Possenti & Gualandris (2002) in- 
vestigated the case of NGC 6752. 



2. The millisecond pulsar in NGC 6397 

NGC 6397 is one of the most promising targets for search- 
ing MSPs in GCs: in fact, it is one of the closest clusters, at 
a distance of 2.6 kpc ± 6% (Reid & Gizis 1998) and prob- 
ably has a collapsed core with hints of mass segregation 
(King, Sosin & Cool 1995). It lies in fourth place in the 
list of GCs ranked according to central luminosity (Harris 
1996). In this cluster only one millisecond pulsar has been 
discovered so far: PSR J1740-5340, a binary MSP with a 
spin period of 3.65 ms (D'Amico et al. 2001a). Its X-ray 
counterpart has been identified by Grindlay et al. (2001) 
among the ~ 20 X-ray sources detected with Chandra 
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Fig. 2. Excess group delays measured in two 128 MHz- 
wide bands centered at 1454 MHz {filled symbols, con- 
nected by a solid line) and at 1326 MHz [open symbols , 
connected by a dotted line) for an event occurring at or- 
bital phases 0.65-0.75 on 2001 March 1 (see Fig. lb). 
The delays in the two subbands are fitted [inserted panel) 
with a straight line of slope 1.21±0.04, corresponding to 
a frequency dependence for the delays (At oc v~'^) with a 
best-fit power-law index 7 = 2.02 ± 0.30 (from D'Amico 
et al. 2001c). 
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2.1. Observed features of PSR .11740-5340 

It displays eclipses at 1.4 GHz for more than 40% of 
the orbital phase. Ten eclipsing systems containing a 
MSP are known (Fruchter et al. 1990; Lyne 1990; Stap- 
pers et al. 1996; Nice, Arzoumanian & Thorsett 2000; 
Camilo et al. 2000; D'Amico et al. 2001a) and in one 
of them (PSR B1744-24A: Nice & Thorsett 1992) the 
eclipses show duration and irregularities similar to those 
of PSR J 1740 -5 340. However this new system is 2 - 3 
times wider (with an orbital separation of ^ 6.5 R0) than 
any other known eclipsing pulsar binary, and has a heavier 
minimum mass for the companion (> 0.19 Mq) than any 
known eclipsing system. 

In addition, the radio signal exhibits strong fluctua- 
tions (delays and intensity variations) over a wide range 
of orbital phases (see Figure 1), indicating that the MSP 
is orbiting within a large envelope of matter released from 
the companion with a high mass loss rate. The frequency- 
dependent behavior of the observed delay events have been 
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studied by splitting the 256-MHz bandwidth available at 
1.4 GHz into two adjacent subbands. Figure 2 shows an 
example of the measured excess delays in the two sub- 
bands. The small panel also shows these measurements 
plotted against each other. The best-fitting straight line 
corresponds to At oc j/-2.02±o.30 j^j, ^j^g excess delays, con- 
sistent with the zy~^ dependency expected if the delays 
arc due to dispersion in an ionized medium. If entirely 
due to dispersion, the excess propagation delays of up to 
~ 3 ms visible in Figure 1 would correspond to electron 
column density variations of ~ 1.5 x 10^* At_3 cm~^, or 
ADM ~ 0.47 At_3 cm~3pc (At_3 being the delay at 1.4 
GHz in ms). This corresponds to a pulse broadening of 

0. 10 P At-3 = 0.36A t-3 ms over the receiver bandwidth 
(to be compared with an intrinsic pulse FWHM of 0.17 P) 
and may be responsible for most of the signal attenuation 
and pulse broadening observed near the delay peaks. 

When combined with the typical duration of the eclipse 
at 1.4 GHz 13 hr). these parameters imply an eclipse 
radius i?ECL 7.5 R©) larger than the orbital separa- 
tion a ~ 6.5 R0 for any value of the orbital inclination 

1. i?ECL is defined as the chord subtended on the orbital 
circumference of radius a by the angle between the orbital 
phase of eclipse ingress (or egress) and the orbital phase 
0.25 (see also Figure 5). 

According to different values of the orbital inclination, 
the values of the Roche lobe radius Rl span the interval 
I.3-2.2R0 (Eggleton 1983). Because i?ECL > a > Rl, 
the matter causing the eclipses escapes the gravitational 
influence of the companion star and therefore must be con- 
tinuously replenished. Assuming isotropic emission, mass 
continuity implies that the companion star releases mass 
at a rate Mc = 47ri?g(~.Lp(i?ECL)w/. Using the values at the 
eclipse radius inferred from the observation of the excess 
delays assuming completely ionized matter (p[i?ECL] ~ 
3 X 10^^^At_3 g cm^^), it results a mass loss rate from 
the companion ~ 1-5 x 10"^"'^ At_3 Vf^s ^Qy~^ where 
Vf^s is the wind velocity at i?ECL in units of 10* cm s~^ 
(typical order of magnitude of the escape velocity from 
the surface of the companion). 

The large position offset of PSR J1740-5340 with re- 
spect to the center of NGC 6397 (r ~ O'.bb) limits the 
unknown contribution to the observed period derivative 
P due to acceleration in the globular cluster potential. 
Thus P can be used to infer the rotational energy loss 
of the pulsar, ij ~ 1.4 x 10^^ ergs s~^. At the distance 
of 6.5 R0, the power impinging on a white dwarf com- 
panion (of typical radius R^d = 0.033 Rq) would be 
Wwd ^ 10^° ergs s~^. This would allow release of an ion- 
ized wind from the degenerate companion at a maximum 
rate M^d.max = W^diRwd/GM^d) ^ 1-6 x lO-^^ Mqj-K 
Thus, for realistic terminal velocity of the wind (w/ ~ 1000 
km s~^) the pulsar flux cannot sustain the observed releas- 
ing of mass, Mc, from a white dwarf companion. 

This conclusion suggested that more likely the com- 
panion could be a bloated main-sequence star or a sub- 
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Fig. 3. Bottom right panel: A portion of a median- 
combined F675W HST WFPC2 image (chip WF4) of 
NGC 6397, centered near the MSP position. The region 
covers about 7" x 7". Star A is the companion to the 
MSP PSR J1740-5340. Remaining three panels: Multi- 
band CMDs for stars detected in a region of the WF4 chip 
(40" X 40") containing the MSP position. The three stars 
found in the vicinity of the MSP error box are marked 
with different symbols (a triangle for star A and squares 
for stars B and C) and labeled with their names in the 
first panel (from Fcrraro ct al. 2001). 

giant, prompting a deep optical search for it. Using a se- 
ries of public Hubble Space Telescope archive exposures 
taken in 1996 March and 1999 April, Ferraro et al. (2001) 
found a bright, ~ 16.6 mag, variable star (hereafter Star 
A) with an anomalous red color and optical variability 
0.2 mag), which nicely correlates to the orbital period of 
the pulsar (see Figures 3 & 4). This variable had been al- 
ready discovered by Taylor et al. (2001); however on the 
basis of the unusual colors and time variability, Taylor 
et al. (2001) had classified this object as a BY Draco- 
nis star. Radius i?c 1-3 — 1.8 R© and effective surface 
temperature Te// ^ 5500 — 5800 K for the companion to 
PSR 1740-5340 can be derived by comparing the color- 
magnitude diagram in Fig. 3 with the isochrones by Sil- 
vestri et al. (1998) and by Vanderberg (2000) for metallic- 
ity [Fe/H]=-2.00 and ages of t = 12- 14 Gyr (compatible 
with the values observed for NGC 6397). 

2.2. The nature of PSR J 1740- 5340 

The red color of PSR-A shows that it is not an undis- 
turbed main sequence star. Its peculiar nature can be un- 
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Fig. 4. Reduced resulting from the fit of two Fourier 
harmonics to the Ha data as a function of the modulation 
period. Small panel: same Ha data folded using the period 
and the reference epoch of the radio ephemeris, and best 
fit (solid line) of the spectral amplitudes of the 1*** and 
2nd harmonics to the data (from Ferraro et al. 2001) 



veiled inspecting the amplitude and the shape of its light 
curves (insert in Figure 4). Two other eclipsing millisecond 
pulsars, both in the Galactic field, have measured light 
curves: PSR B1957+20 (Callanan, van Paradijs & Ren- 
gelink 1995) and PSR J2051-0827 (Stappers et al. 2001). 
In both cases, the optical companions show strong mod- 
ulations, which are readily explained by the heating of 
one side of the companion by radiation from the pulsar. 
The optical maximum is at orbital phase 0.75, where the 
heated side of the companion faces the earth. A similar 
trend, though with a much smaller degree of modulation, 
is seen in 47 Tuc Vopt, the first identified MSP companion 
in a GC (Edmonds, et al. 2001). 

The light curves of Star A are completely different. In 
Figs. 4 & 5, the phase 0.0 is located at the ascending node 
of the MSP orbit; thus at the phase 0.75 the side of the 
companion facing the pulsar is visible. In contrast to the 
other known variable MSP companions, the light curves 
of Star A display there a minimum instead of a maximum 
(see Figure 4). Within the limits in the available orbital 
period coverage, the best-fit light curve of Figure 4 shows 
two maxima and two minima during each binary orbit: 
thus, tidal distorsions appear the more natural responsi- 
ble for this pattern. If this is the correct interpretation, 
the observed modulation can be reproduced only if the 
companion has almost filled up its Roche lobe and the 
orbital plane is nearly edge-on {i ~ 90°). 

In summary, as roughly depicted in Figure 5, this bi- 
nary appears the first example of a system comprising an 
active MSP whose companion keeps on pouring mass in 
the binary system through the inner Lagrangian point of 
its Roche lobe. 



2.3. The origin of PSR J1740~5340 

Three hypotheses have been suggested for explaining the 
origin of this system: 

(a) that Star A is a main sequence (MS) star perturbed 
by the energetic flux emitted from the MSP; 

(b) that PSR J1740-5340 is a new-born MSP, the first one 
observed just after the end of the process of recycling 
(D'Amico ct al. 2001c; Burderi, D'Antona & Burgay 
2002; Ergma & Sarna 2002); 

(c) that PSR J1740-5340 is a re-born MSP, which has al- 
ready experienced many phases of recycling (Grindlay 
et al. 2002). 

Given the relatively large orbital separation, it seems 
energetically difficult to explain an increasing of both the 
effective temperature and of the radius of a main sequence 
companion: however further detailed simulations are re- 
quired before excluding the scenario (a). The hypothesis 
(b) is particularly stimulating because could be the first 
example of a binary in the so-called radio- ejection phase 
(Burderi, D'Antona & Burgay 2002) in which the Roche- 
lobe overfiow stage is not ended yet, but the MSP is al- 
ready active and sweeps away the infalling matter pre- 
venting accretion onto its surface. In particular (Burderi, 
D'Antona & Burgay 2002) Star A could have been orig- 
inally a main sequence star of ~ 1 Mq, whose evolution 
triggered mass transfer towards the compact companion, 
spinning it up to millisecond periods (Alpar, et al. 1982). 
Irregularities in the mass transfer rate from the compan- 
ion, Mc, are common in the evolution of these systems (e.g. 
Tauris & Savonije 1999): even a short decreasing of Mc can 
have easily allowed PSR J1740— 5340 (having a magnetic 
field ~ 8 X 10* G and a rotational period ^ 3.65 ms) to 
become source of relativistic particles and magnetodipole 
emission, whose pressure (i) first swept the environment 
of the "neutron stars (NS), allowing coherent radio emis- 
sion to be switched on (Shvartsman 1970) and (ii) then 
kept on expelling the matter overflowing from the Roche 
lobe of Star A (Ruderman, Shaham & Tavani 1989). For a 
wide enough binary system, once the radio pulsar has been 
switched on, any subsequent restoration of the original Mc 
cannot quench the radio emission (Burderi et al. 2001). 
Remarkably, for the parameters of PSR J1740— 5340, the 
orbital period at which radio- ejection dominates is ~ 39 
hrS; very close to the observed value of ~ 32.5 hrs (Bur- 
deri, D'Antona & Burgay 2002). 

As already noticed, the signal from PSR J1740— 5340 is 
eclipsed for about 40% of the time at 1.4 GHz and strongly 
disturbed along all the orbit (D'Amico et al. 2001c) and 
this suggests a bias against the detection of this kind of 
radio sources. On the other hand, the occurrence of the 
radio- ejection mechanism can strongly shorten the dura- 
tion of the bright X-ray phase in the recycling scenario. 
Hence, if this radio- ejection phase is common in the for- 
mation of the MSPs, the birth-rate of the Low-Mass X- 
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Fig. 5. Sketch of the binary system containing PSR J1740— 5340 at the phase of the egress of the miUisecond pulsar 
from the eclipse at 1.4 GHz. i?ECL represents the eclipse radius, greater than the orbital separation of the system. The 
maximum of the eclipse occurs at phase 0.25 (lower side of the figure). The eclipse ingress is roughly located at phase 
0.05, whereas the eclipse egress is at phase 0.45. 



ray Binaries (depending on their observed number in the 
Galaxy and on the supposed lifetime) should be revised 
upwards. This could help solving the apparent slight dis- 
crepancy between the observed birth-rate of the MSPs and 
that of the Low-Mass X-ray Binaries, which are their sup- 
posed progenitors. 

Finally, hypothesis (c) calls for the high rate of dynam- 
ical interactions between single stars and binaries allowed 
by the crowded environment of the central regions of a 
globular cluster. Due to such interactions, a neutron star 
in a globular cluster could undergo many phases of re- 
cycling, accreting mass from different companion stars ac- 
quired at different times along the cluster evolution. These 
repeated stages of accretion could affect e.g. the evolution 
of the surface magnetic field of the MSPs in globular clus- 
ter and their X-ray luminosity (Grindlay et al. 2002). 

3. The five millisecond pulsars in NGC 6752 

NGC 6752 is one of the globular clusters with the more 
precise distance measurement: 4.1 kpc ± 5%, obtained fit- 
ting the white dwarf (WD) sequence (Renzini et al. 1996). 
It is classified as a core-collapsed cluster with evidence of 
mass segregation (Ferraro et al. 1997). Rubenstein & Bai- 
lyn (1997) estimated a binary fraction for main sequence 
stars in the range 15%— 38% inside the inner core radius 
(< ll")i decreasing to less than 16% beyond that. At least 
6 dim X-ray sources had been identified on the basis of 



ROSAT pointings (Vcrbunt & Johnston 2000), whereas 
a more recent Chandra observation revealed 19 sources 
within the half-mass radius of the cluster (Poolcy et al. 
2002). 

In this cluster a 3.26 ms binary pulsar was first discov- 
ered, and originally labeled PSR J1910-59A (D'Amico 
et al. 2001a). Amplification due to scintillation helped in 
the detection of four additional MSPs in the same clus- 
ter (hereafter PSRs B, C, D, E). All of them are isolated 
with spin periods in the range 4.6 — 9.0 ms. X-ray coun- 
terparts for PSR-D, PSR-C and PSR-B have been de- 
tected (Poolcy et al. 2002, D'Amico et al. 2002) using a 
28700 sec long exposure taken with the ACIS-S detector 
aboard the Chandra X-ray observatory. 

3.1. The three inner miUisecond pulsars 

PSRs B, D and E are located close to the cluster cen- 
ter, as expected as a consequence of mass segregation in 
the cluster (see Figure 6). In particular, PSR— D has the 
third largest period derivative, P = 9.6 x 10~^^, among 
known MSPs (after PSR B1820-30A in NGC 6624 and 
PSR B1821-24 in M28), suggesting that the P is domi- 
nated by the line-of-sight acceleration in the cluster gravi- 
tational field. This interpretation is supported by the large 
negative P values observed for PSR-B and E, which allow 
to directly derive lower limits to the line-of-sight acceler- 
ations, ai/c = P/P = -9.6 ± 0.1 X 10"" s-^ This value 
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Fig. 6. Positions of the 5 millisecond pulsars (arrows) in 
NGC 6752, superimposed on an optical image retrieved 
from the Digital Sky Survey. The circle indicates the half- 
mass radius region {9hm = 115"). 

is the largest known after those of PSRs B2127-I-11A and 
D in the core of M15 (Anderson et al. 1990). Other con- 
tributions to ai/c are all negligible (D'Amico et al. 2002) 
and thus the inferred high values of |a//c| unambigously 
calls for the effects of the potential well of NGC 6752. 

3.2. Central Mass-to-Light ratio in NGC 6752 

A lower limit to the mass-to-light ratio in the inner re- 
gions of NGC 6752 can be derived from the following rule, 
which holds to within ~10% in all plausible cluster models 
(Phinney 1992): 



< 



1 1 G M,_yi{<e±) 



= 5.1 X 10-is^ 



10* Lvopc- 



Here Sy(< 0±) is the mean surface brightness within a 
line of sight subtended by an angle 9± with respect to 
the cluster center, Mcyi{< 0±) is the mass enclosed in the 
cylindrical volume of radius R± = D6± and A4 / Cv is the 
mean projected mass-to-light ratio in the V-band. In Fig- 
ure 7 the curves of maximum \ai/c\ for different values 
of M/Cv have been plotted, using the most recent pub- 
lished brightness profile for this cluster (Lugger, Cohn & 
Grindlay 1995), normalized to the central surface bright- 
ness in the V-band reported by Djorgovski (1993). In par- 
ticular, in Figure 7 the histogram represents the data for 
A^/£v=l-l- the mass-to-light ratio suggested by Pryor 
& Meylan (1993) for the central regions of NGC 6752. 
This greatly underestimates the observed values of P/P. 



The dashed lines are good analytical fits to the observed 
data, scaled according to increasing values of M-/Cy. Only 
M/Cv ^9 can account for the observed \P/P\ of PSRs B 
and E. An even larger M/Cy > 13 is required if PSR-D 
has a negligible intrinsic P/P (see caption of Figure 7). 
It is worthwhile noting that these estimates of the pro- 
jected Ai/C\> are independent of modeling of the cluster 
potential and of distance, excepting the effects of extinc- 
tion. Being E(B-V) very smaU for NGC 6752 (= 0.04 
according to Harris 1996) the latter affects negligibly the 
result. 

It provide the first direct dynamical evidence for a high 
density of unseen remnants in the core of a globular clus- 
ter, suggesting a mass-to- light ratio ^10 in the core region. 
Typical mass-to-light ratios for the central regions of the 
globular cluster span the interval Ai/C\i=l—3.5 (Pryor 
& Meylan 1993). NGC 6256, the GCs at the top of a list 
ordered according to M/Cy (Pryor & Meylan 1993), has 
AA/Cv ~ 6, but this value is quite uncertain due to the 
high level of reddening. Using the same method applied 
to NGC 6752, a much smaller value ~ 3, was obtained by 
Phinney (1993) in the case of M15, a core-collapsed GC 
long suspected to host a central black-hole (see e.g. van 
der Marel & Roeland 1999). 

The nature of the remnants in the core of NGC 6752 is 
not clear as it strongly depends on assumptions about the 
initial mass function of the cluster. However, the possibil- 
ity that they are black holes, or that many stellar rem- 
nants have collapsed into a massive (binary?) black hole, 
is intriguing. 

3.3. The millisecond pulsar in the outskirts of the cluster 

When a coherent timing solution became available, 
PSR— A (the only binary pulsar known in the cluster) 
resulted located at 6 '4 from the globular cluster center 
(R.A.=19h 10™ 51!8, Dec.=-59° 58' 55", eqx 2000, Har- 
ris 1996) equivalent to 3.3 half-mass radii or ~ 57 core 
radii. Hence, it is the farthest objects from the GC center 
known in the catalogue of 41 pulsars in GC with accu- 
rately determined positions. Previously, the largest offset 
of an associated pulsar from a GC center was for PSR 
B2127-hllC (Prince et al. 1991), a member of a double- 
neutron-star eccentric binary in M15, located at ~ 0.8 
half-mass radii. Given this large radial offset, one could 
question the association of PSR— A with the cluster. How- 
ever, the probability of chance superposition of a Galactic 
field MSP with a flux density of S'1400 ~ 0.2 mJy within 
6.'4 of the center of NGC 6752 is at most 10"^ This prob- 
ability is further reduced by noting that all five pulsars 
have very similar DM values. 

The companion to PSR— A is a dwarf (probably de- 
generate) star of mass m* ^ 0.2Mq, moving with period 
-Foib — 0.86 days on a circular orbit around PSR— A (the 
eccentricity is < 10"^, D'Amico at al. 2002). The orbital 
separation is opsr-a 0.0223 AU, and the binary has a 
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Fig. 7. Maximum linc-of-sight acceleration \ai^^^/c\ — 
\P/P\ versus displacement 6 with respect to the center 
of NGC 6752. The histogram represents the prediction 
based on the available optical observations. The dashed 
lines are analytical fits to the optical observations, la- 
beled according to the adopted mass-to-light ratio. The 
filled pentagons represent lower limits to the line-of-sight 
acceleration based on the observation of PSRs B and E. 
The open pentagon shows the value oi P/P for PSR-D as- 
suming a negligible intrinsic P. If various observed scalings 
(Becker & Triimper 1997; Possenti et al. 2002; Grindlay et 
al. 2002) between X-ray luminosity and spin-down power 
for the MSPs are used to estimate the intrinsic P, P/P for 
PSR-D assumes values in the interval between the open 
and the dotted pentagon (from D'Amico et al. 2002). 



gravitational binding energy -EpsR-A ^ lO^^erg. The in- 
ferred magnetic field is ~ 10* Gauss and the characteristic 
age TpsR-A ^ 15 Gyr (given the large offset from the GC 
center, the period derivative of PSR— A is not significantly 
affected by the GC potential well); in short, it looks like 
a "canonical" recycled MSP (Lorimer 2001). 

3.4. The origin of PSR J1911-5958A 

Colpi, Possenti & Gualandris (2002) explored three pos- 
sible explanations for the unprecedented position of 
PSR J1911— 5958A (this is the correct name given to the 
pulsar after the determination of its precise position) : 

(a) that PSR— A, remnant of a massive star, was the com- 
ponent of a truly primordial binary born in the GC 
halo or expelled from the GC core due to the natal 
kick imparted by the supernova explosion; 



(h) that it was ejected from the core via scattering or ex- 
change interaction(s) with cluster stars; 

(c) that it was propelled into the cluster halo by a scat- 
tering event involving a black hole binary (hereafter 
[BH+BH]). 

3.4.1. Case (a) 

Hypothesis (a) is the simplest among the three. It requires 
orbit stability against dynamical friction (having a char- 
acteristic time scale rjf), and a low recoil speed at the 
moment of the neutron star (NS) formation. 

If the binary has been kicked from the globular cluster 
central regions into a highly eccentric orbit (cec ^ 0.9) at 
the time of NS formation, it would remain at the periph- 
ery for a time Tdf ^ 7 x 10^(1. GMq/A/) yr before return- 
ing to the core. Hence, if the supernova took place in the 
central regions, the binary of PSR— A (of mass >2.5Mq 
after the NS has formed) would have been driven back 
to the core by dynamical friction within <1 Gyr. As the 
supernova must have occurred in the first Gyr since GC 
formation, the binary would not have survived in the halo 
until present time. 

For a centrophobic halo orbit (eoc^O.5) the dynamical 
friction timescale Tdf is >10"(1.6Mo/M) yr in NGC 6752, 
implying stability of the orbit of PSR— A and of its progen- 
itor halo binary. However, if PSR— A indeed descends from 
a binary born in the cluster outskirts and following an al- 
most circular orbit in the cluster potential, it must have 
experienced the faintest supernova kick ever observed: in 
fact the escape velocity T4sc is only ~ lOkms"^ in the 
halo of NGC 6752. Conversely, it is known from current 
proper motion measurements of 13 MSPs belonging to the 
Galaxy (Toscano et al. 1998) and from statistical stud- 
ies (Cordes & Chernoff 1997) that MSPs have relatively 
large mean transverse velocities of 85±13kms~^ imply- 
ing peculiar speeds with a 3-D dispersion of 84 kms~^. 
Moreover, theoretical studies estimates that the minimum 
recoil speeds V^cc for a MSP binary is >60kms~^ (Tauris 
& Bailcs 1996). As a consequence, scenario (a) calls for the 
NS of PSR— A being probably formed from the accretion- 
induced collapse of a massive white dwarf (Grindlay et al. 
1988), in which a large kick velocity is not expected. 

A difficulty of this scenario is related with the proba- 
bility that a MSP is formed in the outer halo of a glob- 
ular cluster. It contains only a small fraction of the total 
mass of the cluster and the relatively small stellar den- 
sity does not permit the occurrence of a significant num- 
ber of dynamical interactions. Hence the only available 
channel for the formation of a MSP is just the evolution 
of a primordial binary, similar to the formation scenario 
holding in the Galactic field. Given the estimated number 
-^gai 4 X 10^ of MSPs in the Galactic field (e.g. Lorimer 
2001) and the ratio fj,<5x 10~^ between the mass in all 
the globular clusters to that of the Galaxy, it turns out 
that at most A*"gcs = -/VgaiM/r ^ 20/^ MSPs born from a 
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primordial binary should be hosted by the GCs system 
(comprising ^ 200 objects) at radii larger than r. Here fr 
is the fraction of the total mass in a typical GC which re- 
sides at a distance greater than r from the center. Taking 
r equal to 3 half-mass radii, fr ~ 0.15 and thus A'gcs 3. 
The very low number A'gcs makes improbable the detec- 
tion of a MSP at its offset position, when descending from 
a primordial binary born in the GC halo. 

3.4.2. Case (b) 

Ejection of PSR— A away of the core via a dynamical in- 
teraction with a cluster star calls for two possible scenar- 
ios: (bl) inelastic scattering or (b2) three-body exchange. 
Scattering can occur either before or after recycling as 
it preserves, to a high degree, the circularity seen in our 
binary pulsar. In both cases, a recoil kinetic energy of 
2 — 3 X lO'^^erg is requested for ejecting [NSa+co] (com- 
prising cither PSR— A or the not yet recycled NS and its 
companion) into its halo orbit. This energy must be ex- 
tracted from the [NSa+co] binding energy: inspecting the 
tables of Sigurdsson & Phinney (1993) it is easy to con- 
clude that this energetic requirement imposes an upper 
limit on the orbital period Porb.up of [NSa+co] which is 
<0.4 days. 

At the onset of recycling a "bifurcation" orbital pe- 
riod, Pbit 2 days, separates the formation of a converg- 
ing binary (which evolves with decreasing orbital period 
Porb until the mass-losing star becomes degenerate and 
a short-period binary pulsar is formed) from a diverging 
one (which evolves with increasing Porb forming a long 
period binary pulsar; Tauris & Savonije 1999). Given this 
scenario, PSR— A should result from a converging binary 
and thus its period can not have ever been shorter than its 
current value of 0.86 days and this exclude the scattering 
hypothesis (bl). 

Exchange (scenario (b2)) alters the nature of the inter- 
acting system and makes the new binary rather eccentric 
and with a longer orbital period. Thus due to the very long 
circularization time between compact stars, a three-body 
exchange must occur before recycling (reminding that the 
upper limit on the eccentricity of the orbit of PSR— A is 
10~^). The simplest and most probable case in this sce- 
nario is that of a main sequence binary [MS-|-ms] (MS 
heavier than ms) interacting with the NSa (slowly rotat- 
ing not yet recycled progenitor of PSR— A) flying by (given 
the similar masses of MS and NSa, the results are almost 
unaffected when swapping their roles). In such a three- 
body encounter the energetic requirement on the recoil 
kinetic energy of [NSa+co] can be accomplished, provided 
that the original main sequence binary is tight enough. On 
the other hand the tables of Sigurdsson & Phinney (1993) 
show that after the exchange, the [NSa+co] binary has an 
orbital period suspiciously close to the currently observed 
value, considering that recycling (and hence shrinking of 
the orbit) has not occurred yet. Even more disturbing is 



the value of the characteristic time (~ 18 Gyr > age of 
the cluster) for an interaction of this kind to take place in 
the core of NGC 6752. 

3.4.3. Case (c) 

Partly motivated by the unusual high mass-to-light ratio 
discovered in the inner region of NGC 6752 (sec previous 
section), Colpi, Possenti & Gualandris (2002) investigated 
a more exotic scenario, involving the presence of a binary 
black-hole in the centre of the cluster. In effect, dynami- 
cal studies (Portegeis Zwartz & McMillan 2000) show that 
some clusters could still retain one binary, the last which 
has not been kicked out. This would be the consequence of 
repeated exchange interactions and ejections (Sigurdsson 
& Hernquist 1993; Portegeis Zwartz & McMillan 2000), 
occurring early in the cluster lifetime among the initially 
rich population (of few hundreds) of BHs, born from the 
most massive stars in the cluster. As they sink toward 
the center by dynamical friction (on a timescale <10^ yr), 
BH— [BH+BH] interactions become overwhelmingly im- 
portant. Binaries hardens progressively in virtue of these 
encounters and eventually leave the cluster due to recoil. 
As a result, many clusters lose their original BHs on a 
time scale of ~ Gyr. But the outcome of the last encounter 
(which could involve two single BHs or a single BH and 
a binary BH system) depends on the specific and unpre- 
dictable parameters of the interacting corpses and thus 
could either eject all the 3 bodies or leave a single BH or 
leave a binary BH. Following Miller & Hamilton (2002a), 
the supposed binary that avoided ejection should likely 
have unequal masses and could harbor the heaviest BH 
that NGC 6752 ever had. 

The scenario (c) is fully consistent and viable provided 
that: 

i the [BH+BH] binary that survived the early ejection 
from the GC, later avoided also coalescence by emis- 
sion of gravitational waves (GWs), at least for a time 
row > 10 Gyr; 

ii the separation of the [BH+BH] binary is sufficiently 
large to have prevented ejection from the cluster due to 
a single interaction with a much lighter stellar intruder; 

Hi the hardening of the [BH+BH] binary due to the en- 
counters with the bath of the cluster stars impinging 
on it does not reduce its orbital separation enough for 
triggering coalescence by emission of GWs (point i) or 
ejection (point ii); 

iv during the scattering interaction, the [NSa+co] binary 
is not ionized by the [BH+BH] binary tidal field; 

V the [BH+BH] binary can have transferred enough ki- 
netic energy to propel PSR— A into the cluster halo; 

vi both the probability of occurrence of the scenario (c) 
and the probability of detection of PSR— A at the 
present offset position arc not negligible. 
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Fig. 8. [BH+BH] binary separation a against mBn/Af© 
for A/bh = SOAf©; acw (heavy solid line), Orec (heavy 
dot-dashed lines) and Ogcat (heavy dashed line) are com- 
puted from the related equations presented in the para- 
graph 3.4.3. In particular arcc, is computed for the case of 
light m = IOA/q black-hole intruders and for the case of 
intruders as massive as the secondary black-hole, namely 
m = mBH- The dotted line gives the tidal disruption ra- 
dius for [NSa+co]. The shaded area indicates the values 
of a for ejection of PSR— A during the last Gyr, compati- 
ble with its detection (from Colpi, Possenti & Gualandris 
2002). 

Indicating with A/bh and ttibh the masses of the two 
BHs in the binary, with AfBH.t and /ibh their total and 
reduced mass, and with m the mass of an intruder, a 
[BH+BH] binary (with typical eccentricity of 0.7) satis- 
fies the condition i) ii its separation a exceeds a critical 
value ogW) that is 



a > acw 0.4 



A'/, 



BH,t 



1/2 



A^BH '''GW 



1/4 



AU 



IOOMqJ \10Mq m^'^yr , 
On the contrary the condition ii) imposes the inequality 



a > flrcc ^ 0.01 I : 



lOOAf© 35kms" 



lOAf© A 1OA//0 AfBH.t K,e 



AU . 



For the parameters of NGC 6752, the (very approximate) 
number of stars impinging onto the target binary over a 
time t - 10^° 



yr IS 



AC < 2 • 10^ 



whereas the minimum required number of encounters for 
triggering the coalescence of the [BH+BH] binary due to 



gravitational waves emission is 

A^GW 5 • 10^ . 

Hence condition Hi) is accomplished. The condition iv) 
requires that the distance of closest approach between the 
two binaries satisfies the relation 



a > rtid ^ 0.2 



/ Af] 



BH,t 



1.6A'/, 







\ IOOA/q mpsR-A 



1/3 

AU . 



The condition wj on the conversion of the binding energy 
of the BH binary into kinetic energy for the ejection (with 
initial velocity Vcj) of the [NSa+co] binary can be accom- 
plished if 



a < a.. 



BH 



L3 



/ Mbh 
\10Mp 



35 km s 



AU 



Finally, one can estimate a rate for MSP ejection TZmsp of 
~ 0.5 Gyr~^, implying a probability of detection of about 
~ 50% considering a life-time for the current orbit of 1 
Gyr. The value of TZmsp can be derived on the basis of 
the flux of stars on the [BH+BH] over the cluster lifetime 
(see the equation for Af*), assuming that about 1/30 of 
these encounters should have involved a NS and that ~ 5% 
of the NSs is in the state of a MSP. These figures are 
calculated adopting a Salpeter-like initial mass function 
(giving ~ 1% of NSs) and a mass segregation enhancement 
of ~ 6 in the density (Sigurdsson & Phinney 1995). If 
the IMF is fiatter (as suggested by the observations of 
D'Amico et al. 2002), this rate could be higher, enhancing 
the likelihood of this scenario. Once ejected into an halo 
orbit, [NSa+co] spends only ~ 1/50 of its fife in the dense 
GC core and thus has negligible probability to undergo 
further interactions. 

Figure 8 summarizes the results of all the previous con- 
siderations, showing that a dynamical interaction with a 
[BH+BH] binary in a mass range of ~ 3 — lOOAf© and sep- 
aration of ~ 0.3-2 AU can propel PSR— A into the cluster 
halo, being consistent with the conservation both of en- 
ergy and linear momentum. In particular an unequal mass 
binary with a heavy >lOOAf0 BH and a lighter compan- 
ion ~ IOMq is preferred, as it can form in a GC (Miller & 
Hamilton 2002b) and can provide the right statistics for 
explaining PSR— A. 

4. Conclusion 

During the last two years, the detection of 15 millisec- 
ond pulsars in six globular clusters which contained no 
previously known pulsars (Possenti et al. 2001; Jacoby et 
al. 2002) broke a 7 years long hiatus in such discoveries. 
More recently, the first stage of timing observations for 
most of the new objects has been completed, opening the 
possibility of investigating both their properties and the 
characteristics of the globular clusters in which they are 
embedded. Here I have presented the first results of these 
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studies in the cases of the globular clusters NGC 6367 
and NGC 6752. They confirm that the efforts devoted to 
the search and timing of millisecond pulsars in globular 
clusters, when combined with observation at other wave- 
lengths, can deserve a significant and often surprising sci- 
entific pay-back. 
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